The states l 1 ^, l 1 !!^ and 2 1 £+ of Sr 2 studied by Fourier-transform spectroscopy 



O 

> 
O 

in 

(N 

S 1 

i 

S 

<D 

O 
c/5 
O 

CZ5 
>> 
43 



> 

m 

On 
O 



Alexander Stein, 1 Horst Knockel, 1 and Eberhard Tiemann 1 

1 Institut fur Quantenoptik, Leibniz Universitat Hannover, Welfengarten 1, D-30167 Hannover, Germany 

(Dated: January 20, 2013) 

A high resolution study of the electronic states 1 E+ and 1 n n which belong to the asymptote 
4 X D + 5 1 S and of the state 2(A) 'St , which correlates to the asymptote 5 1 ? + 5 1 S, is performed 
by Fourier-transform spectroscopy of fluorescence progressions induced by single frequency laser 
excitation. Precise descriptions of the potentials up to 2000 cm -1 above the bottom are derived and 
compared to currently available ab initio calculations. Especially for the state 1 X E+ large deviations 
are found. Rather weak and local perturbations are observed for the states l 1 ^ and 2 E„, while 
a strong coupling of the state 1 1 E^ to the component Q — 0^ of the state l 3 n«, which belongs to 
the asymptote 5 3 P + 5 1 S, is indicated. 

PACS numbers: 34.20. Cf, 31.50.Df, 33.20.Kf, 33.20.Vq 



I. INTRODUCTION 

Our current spectroscopic work on the Sr2 molecule is 
motivated by the interest in ultracold Sr 2 dHU, where 
the desired information on the molecule for planning ex- 
periments was obtained up to now from ab initio calcula- 
tions [J]. In our earlier work on Ca 2 [5J], an isoelectronic 
molecule to Sr2, we observed that the available ab initio 
calculations are not as precise as they are e.g. for alkali 
molecules. Thus we want to test the accuracy of the the- 
oretical calculations which are available for Sr 2 0, H, 0K 

In our previous publications on this molecule @, 0] 
we mainly concentrated our effort on the ground state 
X^^E^ . For the present work we performed extensive new 
measurements to increase the knowledge on the excited 
states. Figure [1] gives an overview of the electronic states 
which arc important in the current work. The potentials 
are taken from the ab initio work [|| which shows 
the closest approach to our present observations. 

We report on the successful observation of the state 
which is presumably only known from the spec- 
troscopy in rare gas matrices [Tj] - [l3j , on the observation 
of the state l 1 !^, which, to our knowledge, was not ex- 
perimentally known before, and on an extended data set 
of the state 2 1 £+, which was in the past the mostly in- 
vestigated excited state, called A state [3 EH, and was 
already used for exciting the fluorescence to the ground 
state in our first paper on this molecule 

The paper is organized as follows: Section [TT1 describes 
the experimental methods, Sec. IIIII describes the ob- 
tained data set and discusses the measurement uncertain- 
ties, Sec. [IV] presents the resulting sets of Dunham co- 
efficients and potential descriptions and compares them 
to the currently available theoretical calculations, while 
Sec. IVl gives conclusions and a short outlook. 



II. EXPERIMENT 

The experimental procedure and the apparatus used 
are the same as in [8f. Strontium is filled into a stainless 
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FIG. 1: (Color online) Overview of the electronic states 
needed for interpretation in this paper. The potentials stem 
from the ab initio calculations [tj. [lol]. 



steel heatpipe, which is heated to a temperature of about 
1220 K with 20 mbar of argon as a buffer gas. To keep 
the optical path free from condensing Strontium which 
grows as crystals at both ends of the heatpipe, the oven 
has to be moved every two hours to melt and bring the 
condensed Sr back to the heated zone. This means, that 
the heated cell is not working perfectly as a heatpipe 
where the condensed substance flows back to the heated 
region. 

A Fourier-transform spectrometer of type Bruker IFS 
120 HR resolves fluorescence progressions induced by a 
cw laser beam. For most spectra a resolution of 0.05 
cm -1 was chosen and typically about 20 scans of the 
spectrometer were averaged. This method is most effi- 
cient if fluorescence to the wanted state can be observed 
by the excitation of higher lying states, because each flu- 
orescence progression contains typically information on 
a larger amount of different rovibrational levels (with 
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AJ±1 or 0) of the lower state. However, for the ex- 
cited state only one excited rovibrational level is involved. 
Fluorescence to states other than the ground state was 
searched without success. Even with a high power Ar + 
laser (2.5 W), which was used to investigate the existence 
of possible very weak transitions and two-photon excita- 
tions, no fluorescence to states other than the ground 
state was observed. This approach was also tried unsuc- 
cessfully with UV lines of that laser. 

An alternative for a more efficient acquisition of a large 
amount of data could have been absorption methods. It 
was tried to record the absorption spectrum of Sr 2 with 
a white light source, but even with long averaging times 
(1000 scans which takes more than 30 hours) only signal- 
to-noise ratios of up to 5 could be reached. Moreover, 
the comparison with simulated spectra based on the po- 
tential regions already known from our earlier work Q 
showed that the visible absorption peaks do not stem 
from single strong lines but from an overlap of many 
lines. This finding is in good agreement with the recorded 
fluorescence spectra, where the laser, independent of its 
exact excitation frequency chosen within the frequency 
range of a strong electronic transition, excites each time 
more than a dozen of rovibrational levels. The reason 
for this is due to the large mass of the Sr dimer together 
with the high working temperatures implying that almost 
all existing rovibrational levels (about 5400 for the ref- 
erence isotopologue 88 Sr 2 ) of the relatively flat ground 
state potential (the binding energy is 1081.64(2) cm -1 
@) are significantly thermally populated. Other meth- 
ods as polarization labe ling spectroscopy were already 
tried by the authors of |15l| w ithout success because of 
too weak molecular signals [16j . 

Due to the high amount of simultaneous excitations 
the method of recording fluorescence induced by the di- 
rect laser excitation of the wanted electronic states is 
not as inefficient as it seems to be at first glance. The 
recording of typical 20 scans with a resolution of 0.05 
cm -1 takes about 20 minutes yielding information on all 
simultaneously excited levels. In addition, the strongest 
progressions sometimes are accompanied by a few colli- 
sionally induced rotational satellites. During the record- 
ing process the experimentator can already use a widely 
automated software to assign the lines of the previously 
recorded spectrum and incorporade the obtained transi- 
tion frequencies and level energies to improve the predic- 
tions for the next excitation frequencies. 



A. Observation of the state 1 E„ 

With a titanium sapphire laser in the range from 12000 
cm -1 to 14000 cm -1 we collected a sufficient amount of 
data for a precise description of the bottom of the po- 
tential of the state l 1 !^ . Only v' = for this state was 
not observed because of too small Franck-Condon overlap 
with the vibrational manifold of the ground state. During 
this investigation we also tried a dye laser operated with 



DCM dye, searching the frequency range from 15250 to 
16700 cm -1 for higher excited vibrational levels of this 
state. Nine progressions were observed which were ex- 
cited close to 15250 cm -1 and originate from three neigh- 
boring vibrational states of that electronic state. These 
levels with term energies around 15500 cm -1 above the 
potential minimum of the ground state (but below the 
atom pair asymptote 1 So+ 3 Pi at 15586 cm -1 ) show sig- 
nificant local perturbations which became obvious from 
plots of these term energies as a function of J (J + 1). 
Additionally, 18 progressions were observed with excita- 
tion frequencies around 16250 cm -1 belonging to seven 
different vibrational levels. They form a structure of frag- 
mentary bands in the term energy region between 16790 
cm -1 and 16890 cm -1 . The line intensities in the center 
of these bands are as strong as the strongest ones ob- 
served for low vibrational levels of 1X1^. The observa- 
tion of fragmentary bands clearly indicate predissociation 
to the 3 P asymptotes between 15586 cm -1 ( 1 S + 3 Pi) 
and 15980 cm -1 ( 1 S + 3 P2)- Though the frequency in- 
terval from 16000 cm -1 to 18500 crn^ 1 was intensively 
investigated using dye lasers operated with DCM, Rho- 
damine 6G and Rhodamine 110, only three progressions 
belonging to levels outside the mentioned energy interval 
were observed and excited with laser frequencies close to 
18250 cm . These three levels do not fit into the now 
well known potential energy curves of the other two elec- 
tronic states in this frequency region (1 n u and 2 1 £+) 
and the progressions have line intensities which fit to 
Franck-Condon factors expected for highly excited lev- 
els of the state l 1 !^. Since also the observed rotational 
and vibrational spacings fit into the expectations for this 
state, it seems to be obvious that the levels belong to a 
second region with fragmentary bands of this state. 



B. Observation of the state l 1 !^ 

During the investigation of the fragmentary bands of 
the state l 1 ^ one progression of Q lines was observed 
which was later assigned tov' — l of the state l 1 !^ corre- 
lating to the asymptote 1 So+ 1 D2, see Fig. [T] This state 
was further investigated, first by searching for neighbor- 
ing rotational levels J' of the same v' and then by search- 
ing for neighboring v 1 levels. For this search theoretical 
B e and oj e from [6j were used until enough levels were 
available to fit a preliminary set of Dunham coefficients 
which then were used for further extrapolations. The ef- 
fort in time for the investigation of the state l 1 !!,, was 
significantly higher than for the states because the 
observed line intensities of the Q progressions belonging 
to this state are typically a factor of seven lower than the 
line intensities of the P, R progressions from the states 
Because of this fact the chance of finding acciden- 
tally excited transitions is low and the precision needed 
for successful predictions of new transitions is high. Vi- 
brational levels of the state l 1 !^ were investigated up to 
v' = 22 using dye lasers with DCM, Rhodamine 6G and 
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Rhodamine 110 reaching excitation energies up to 18500 
cm . 



C. Observation of the state 2(A) 1 E+ 

Since the transition strengths from the state 2(A) 1 E+ 
showed to have on average a similar amplitude as for 
the state l 1 ^, it was also possible to collect a sufficient 
amount of data for a precise description of the potential 
up to 2000 cm -1 above its bottom using dye lasers with 
Rhodamine 6G and Rhodamine 110. Nevertheless, the 
overwhelming amount of data obtained for this state is 
a side product of the investigation of the ground state 
0, [9| and the intensive search for energy levels of the 
state l 1 !^ as mentioned above. Using the visible lines 
of an Ar + laser a large amount of time was investigated 
to search for highly excited levels of this state, but at 
the end this was not very successful. A few progressions 
were found by exciting with the 514 nm and the 496 nm 
lines. To identify few weak progressions excited with the 
488 nm line much longer averaging time and the use of 
a single-mode assembly for the Ar + laser were necessary. 
For the more blue lines of the laser no fluorescence pro- 
gressions were found. This strong decrease of the fluores- 
cence intensities towards higher excitation energies can- 
not be explained by low populations of those ground state 
levels with best Franck-Condon overlap since the ground 
state population changes only by a factor of two from 
the potential minimum to the asymptote for the high 
working temperature of 1220 K. Because of the high level 
density in the long range regions of both the excited and 
the ground state the number of simultaneous excitations 
strongly increases for high v 1 and because the maximum 
total absorption of the Sr2 molecules in the heatpipe is 
only of about 10% the strength of the total fluorescence 
should strongly increase. But we observe that it stays 
roughly constant and the individual progressions become 
too weak for a reliable analysis. The most likely explana- 
tion is again predissociation which could be caused by a 
growing Franck-Condon overlap of the higher vibrational 
states v' with the continuum levels of the lower lying 
triplet and singlet states. 



III. DATA SET 

The data helds of the three excited states are given in 
ngurc H and table Q] gives the numbers of levels for the 
different electronic states and isotopologues. The iso- 
topologue 86 Sr2 is only observed through levels with low 
quantum numbers v and J where it cannot be distin- 
guished from the isotopologue 84 Sr 88 Sr because of almost 
equal reduced masses. There are no independent obser- 
vations of the isotopologue 84 Sr 88 Sr through levels with 
even J', which do not exist for states of type of 
the homonuclear isotopologue 86 Sr2. Both isotope com- 
binations have also similar natural abundances, thus we 
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FIG. 2: (Color online) Overview of the observed energy levels 
of state (a), state l 1 n u __^b) and state 2 X E^ (c), including 

all measurements of ref. [ 



TABLE I: Number of observed levels for the investigated elec- 
tronic states for various Sr2 isotopologues. The last column 
gives the number of ground state levels observed by fluores- 
cence from the excited states. 





abundance 


1±E+ l 1 ^ 2 X E+ 


X*E+ 


84 Sr 88 Sr 


0.92% 


j. 2 


70 




86 Sr 2 


0.97% 






86 Sr 87 Sr 


1.38% 


2 


42 


87 Sr 2 


0.49% 


blended by 86 Sr 


88 Sr 


86 Sr 88 gr 


16.28% 


37 3 223 


4485 


87 Sr 88 Sr 


11.56% 


23 1 142 


2955 


88 Sr 2 


68.19% 


272 149 675 


4890 


total 


99.81% 


332 153 1044 


12442 
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expect that the observed progressions are the sum of con- 
tributions from both isotopologues which are excited by 
the same laser frequency. By the same reason the pro- 
gressions from the isotopologue 87 Sr2 for low quantum 
numbers should always be blended by the stronger (fac- 
tor of 33) progressions from the isotopologue 86 Sr 88 Sr. 
Thus it was never observed separately during these mea- 
surements. Of the 153 observed levels of the state l 1 n u 
only 13 are e levels, all remaining rovibrational levels are 
f levels. The data set for the state 2 1 E+ published in [H 
consisted only of 260 levels compared to the 1044 used for 
the current analysis. The levels of all individual states 
are derived from altogether 60000 assigned lines of which 
15700 were taken from our previous works [1, H|. 

The assignment was done using a specially developed 
software as described in [aQ (in a prior version in [l7|). 
The assignment function mostly used for the current 
work needs as input one selected line from the loaded 
spectrum which is assumed to belong to a progression and 
the excitation frequency, where it takes as default value 
the frequency of the strongest line (if the user does not 
set another frequency) . The program then runs through 
all v" and J" assignments for all possible isotopologues 
of all molecules of which the data of the rovibrational 
manifolds are loaded (for this work only data for Sr2 was 
loaded). The first test done for each possibility of as- 
signment is if the laser frequency fits into the progression 
corresponding to the quantum numbers, thus if the as- 
sumed excited level could have been excited by the laser 
frequency. If this is the case the spectrum is searched for 
additional lines which could belong to the progression. 
The final rating criterion for Q progressions is the high- 
est number of additionally found 'good' lines, i.e. lines 
which do not overlap with other lines and have a signal- 
to-noise ratio better than three. In the case of a doublet 
progression the program selects the progression with the 
highest number of 'good doublets', i.e. doublets consist- 
ing of two good lines with intensity deviations from each 
other of less than 30%. This selection is presented to the 
user, who then has to decide if he wants to save or dis- 
card the suggested assignment. This assignment function 
is a member of a higher order subroutine which automat- 
ically applies the function to all yet unassigned lines in 
the current run of evaluation. This macro routine was of- 
ten run during the nights between measurement days to 
search for progressions overlooked by the user in the less 
automated assignments during the measurements. But 
all automatic assignment are checked in the potential fit 
later, which gives the stringent consistency proof of the 
whole dataset. 

The absolute vibrational assignment of the states l x £„ 
and l 1 !^ was first derived by the comparison of the ob- 
served relative line intensities to Franck- Condon factors 
calculated from preliminary RKR potentials obtained for 
different vibrational assignments. A vibrational assign- 
ment by counting the minima of the intensity envelopes of 
the observed progressions was not possible for these two 
states because the positions of the potential minima are 



too strongly shifted towards small internuclear distances 
compared to the ground state and thus not all minima 
of the intensity envelopes are visible. The vibrational as- 
signments are later confirmed by the observation of fluo- 
rescence progressions from the weaker isotopologues with 
the correct frequencies predicted from mass scaling. 

The vibrational assignment of the state 2 1 £+ was 
taken from our previous work Q in accord to the re- 
sult in [l5| and is later confirmed here by the methods 
described above. 

For all measured lines we assume an uncertainty of 0.01 
cm -1 for absolute line frequencies which equals the un- 
certainty of the Fourier spectrometer specified for abso- 
lute frequencies (in contrast to the specified uncertainty 
for frequency differences of 0.001 cm -1 ). The Doppler 
half width for Sr2 at a temperature of 1220 K is 0.02 
cm , but because of the weakness of the observed Sr2 
fluorescence it is less probable to detect progressions ex- 
cited with laser frequencies in the wings of the line pro- 
files. So the uncertainty estimation of 0.01 cm -1 should 
be justified. 

IV. RESULTS 

As already mentioned in the above sections, to each 
level of the excited states belongs a fluorescence progres- 
sion with individual rovibrational levels of the ground 
state which are employed for the assignment. The en- 
ergy of the excited level is then calculated by averaging 
the sums of the line frequencies and the energies of the 
corresponding ground state levels. Thus, the level energy 
and the quantum number J' due to the selection rules are 
fixed by the very precise ground state potential [9|. The 
lists of the observed excited levels can be found in the 
additional online material. 

The rovibrational level scheme of an electronic state 
with electronic quantum number can be described by 
the conventional Dunham expansion 18] in many cases: 

Kj - E y» (yf ) Wk (v + 1) 1 ■ W + 1) - k • 

The Yik are the Dunham parameters of the chosen ref- 
erence isotopologue 88 Sr2 with reduced mass ^o, while 
/ij is the reduced mass of the isotopologue i, for which 
the energy ladder is considered. This simple approach is 
used for a first consistency check and for the convenience 
of predicting new transitions. Later the information will 
be concentrated in the derived potential energy curve of 
the electronic states. 

A. The state 1 X S+ 

The reduction of all observed transitions leads to 332 
different levels of state l^Szt. 
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TABLE II: Dunham coefficients for the state l 1 !)^, for the reference isotopologue 88 Sr 2 with 1 < v' < 10 and J' < 220. All 
values in cm" 1 . The energy scale is referred to the potential minimum of the ground state. 
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Table [TT| gives a set of Dunham coefficients derived for 
the quantum number range 1 < v' < 10, J' < 220 of the 
state Attempting to produce a set of coefficients 

for the full range of observation it soon became obvious 
that for higher v' and J' the levels are strongly perturbed. 
The value given for Yqq is the sum of the original Dun- 
ham correction Y 00 [18[ and the electronic term energy T e 
calculated with respect to the potential minimum of the 
ground state X 1 !^, which is 1081.64(2) cm" 1 below the 
asymptote 1 So+ 1 So- The coefficient set reproduces the 
196 levels within the given quantum number range with 
a weighted standard deviation of a = 0.59 which might 
indicate that the assumed measurement uncertainty of 
mostly 0.01 cm" 1 might be too large. The validity of the 
Kratzer relation (deriving the centrifugal distortion (Y02) 
from the rotational (loi) and vibrational constant (Vio): 
Y 02 ~ — $Yqi/ Y\q), which holds within 2%, indicates that 
the minimum region of the potential is not very much af- 
fected by the perturbations. The number of digits given 
in table [TTI is adjusted to allow a precise reproduction of 
the observed energy levels, i.e. better than 0.01 cm" 1 . 

Since the amount of currently available data in the 
perturbed potential region is not sufficient for a deper- 
turbation analysis we decided to try a single potential 
fit. Because of the expected appearance of an avoided 
crossing (see figure [IJ a potential fit using an analytic 
potential representation as in Q will not give satisfying 
results when all observed levels up to v' — 20 were in- 
cluded. Thus we used a more flexible potential descrip- 
tion with cubic splines based on the algorithm of [19| 
(also described in [2(|) to fit a potential, which clearly 
shows the turn over to an avoided crossing reaching the 
energy where the ab initio calculations [J @, 0] predict 
the state l 3 n„. Obviously the adiabatic picture will be- 
come more appropriate for the current situation. Ac- 
cording to this picture the derived potential energy curve 
describes the state 1 0+ , see figure [3J (full lines) , which 
coincides in the minimum with the state l 1 ^]^", but dis- 
sociates to the atomic asymptote 3 Pi + 1 So- The central 
part of the potential (Ri < R < R a ) is described by sam- 
ple points interpolated with cubic spline functions (see 
e.g. [2lj]) and is connected at the end points Ri and R a 
to short and long range extensions with continuously dif- 
ferentiable connections to the extension functions (con- 
trary to [19, 20j). The short range extension at the inner 



repulsive wall (R < Ri) has the form 



Vi(R)=A + 



B 

R™ ' 



and the long range part (R > R a ) is given by 



(2) 



(3) 



where Uoo is the asymptotic energy of the pair Sr Pi + 
Sr 'I'.;,. 

It is important to note that the coefficients Cq and Cg 
are only used to obtain a reasonable potential shape in 
the region where no or only few data are available and 
they should not be interpreted as long range coefficients, 
which would allow conclusions on atomic properties. 

Figure [3] shows the derived potential curve as thick full 
line. The position of the avoided crossing in the resulting 
potential, around 4.5 A, clearly indicates that the per- 
turbing state is the state l 3 n 0ll . The potential describes 
the data set (measurement precision of 0.01 cm -1 ) with 
an accuracy better than a few tenth of cm" 1 in most 
cases. The reason for such large residuals should be re- 
lated to perturbations by the spin-orbit coupled state 
for which at least one level is observed (likely v' = of 
£1 = 0) and by the £1 — 1 component of the state 1 3 II U 
with its rotational coupling to 3 no, t . For the potential fit 
the data belonging to higher lying adiabatic potentials of 
the coupled system marked in figure[3jwith 'low v", 'weak 
progressions' and 'fragmentary bands' are currently not 
used but may be helpful for a future deperturbation anal- 
ysis. 

The potential coefficients are given in table IIII1 For 
obtaining level energies within the lowest ten vibrational 
states it is much better to use the Dunham coefficients 
given in table [Til or an alternative potential function con- 
tained in the online material, which describes the same 
reduced data set with v' < 10 as the Dunham coefficients 
do with a weighted standard deviation of a — 0.65. 

The potential curve is certainly precise enough for 
a comparison with ab initio potentials being the main 
goal of the present work. Figure [3] compares graphically 
the resulting potential to the ab initio calculations from 
Boutassetta et. al. @, (dashed line), Czuchaj et.al. 
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TABLE III: Potential coefficients for the state 1 the po- 
tential energies are calculated with respect to the potential 
minimum of the ground state. The given long range coeffi- 
cients are only derived for a proper extrapolation to the dis- 
sociation asymptote 1 So+ 3 Pi- 



R [A] Energy [cm x ] 



3.464536 
3.546908 
3.603682 
3.635081 
3.702534 
3.766714 
3.821524 
3.886870 
4.001479 
4.061114 
4.147573 



14161.788 
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13407.698 
13283.387 
13068.393 
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12843.185 
12792.860 
12821.585 
12867.134 
12963.159 



R< Ri = 3.464536 



A 10653.785 cm" 1 
B 6066396 cm -1 A 6 
n 6 



R[A] 
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4.209351 
4.259994 
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4.446272 
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4.542847 
4.622617 
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4.968028 
5.119717 
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13728.271 
13873.598 
14042.967 
14256.570 
14483.744 
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FIG. 3: (Color online) Comparison of the potential of the 
state l J Ej with the ab initio curves from Boutassetta et.al. 
[I GH, Czuchaj et.al. and from Kotochigova 4] calculated 
for the states 1 X E^ and 1 3 II U . 



.7] (dotted) and from Kotochigova [I[ (dash-dotted) de- 
rived for the states and 1 3 II U . The potential data 
from 0, [l"o| were originally given in hartree as energy 
unit and thus converted to cm" 1 for this comparison 
(1 hartree = 219474.63 cm" 1 ). We adjusted the en- 
ergy offset of the complete potential set from [|| [Io| such 
that their ground state asymptote lies at the dissociation 
asymptote as given in @. Usually it is more appropri- 
ate to take the potential minimum as the reference point 
for comparisons to ab initio calculations, but in the cur- 



rent case forQ no ground state potential is available and 
for Q and [J] the deviations in the dissociation energies 
of the ground state potentials are negligible small com- 
pared to the excited states. The potential curves from 
Q and [|| were originally referenced to the ground state 
asymptote and thus the energy offsets for these potentials 
were shifted by the experimental ground state dissocia- 
tion energy reported in [9( . In figure [3] one can clearly 
see that all the potentials calculated for the state 1 3 IT U 
agree quite well (up to an internuclear distance of 4.5 A) 
while the potentials calculated for the state l 1 ^ show 
deviations as if they were different electronic states. It 
might be surprising, that the earliest work by Boutas- 
setta et.al. [6 1 comes closest to our experimental results 
and overestimates the binding energy of this state by 
about 400 cm" 1 , while the two recent works by Czuchaj 
et.al. and Kotochigova [3| underestimate it by about 
3000 cm" 1 . For the strong deviation of the result from 
we note, that their potential looks like the 1 1 A. U 

potential from Q and their l 1 A ll potential looks like the 
l 1 !]^ potential from [f|. If they were wrongly assigned 
they would only have deviations of the small magnitude 
as for Kotochigova Q did 'full-relativistic' calcula- 
tions, published adiabatic potential energy curves and 
wrote that 'various avoided crossings' between the low 
0+ states exist. Thus it is not completely clear which 
of the reported potential curves should be used for this 
comparison. Because the published pictures suggest that 
all avoided crossings are at the inner repulsive poten- 
tial walls, while the states look well separated and well 
shaped in the regions from the asymptotes to the poten- 
tial minima, we took the state 3 0+, which is the state 
with = from the asymptote X D + 1 S. In any case Ko- 
tochigova did not predict the observed avoided crossing 
between the states l 1 !!^ and l 3 rio u and thus definitively 
the result in Q shows a significant conflict with the ob- 
servations. 

One could ask if our assignment of the electronic state 
is correct and if we really observed the state l 1 ^. We 
observed strong line intensities and only transitions with 
A J = ±1 and no Q lines, thus it is extremely unlikely 
that we observed anything different from a state of type 
1 E+. The £1 = 0+ component of a triplet state would be 
accessible because of the strong spin-orbit coupling of the 
Strontium atom but should have much lower transition 
dipole moments to the ground state than a 1 E+ state. 
Additionally, our present observations are very consistent 
with the experimental observations for the quite similar 
molecule Ca 2 @. 

Contrary to the calculations for the state l 1 ^^ t ne a ^ 
initio potentials [1, [f| 0] of the state l 3 n u will fit to the 
observed position of the avoided crossing between these 
two states. Only the long-range behavior of the result 
from [H would correspond to an unrealistic large C3 co- 
efficient which is in contradiction to the measurement 
results from ultracold ensembles (e.g. @)- 

Table IIVI summarizes for the state l 1 ^ the exper- 
imentally determined spectroscopic constants T e , equi- 
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TABLE IV: Comparison of spectroscopic constants of the 
state FE^ and the isotopologue 88 Sr2. The energy T e is cal- 
culated with respect to the potential minimum of the ground 
state, our values for T e and R e are taken from the potential 
fit, while the values for u e and B e are taken from the Dunham 
fit. 



Source T e [cm 1 


Re [A] UJ e [cm V 


B e [cm- 1 ] 


This work 12796(2) 


3.95(1) 80.71(3) 


0.024794(2) 


[6] 12363 


3.850 79 


0.0259 


[7] 15792 


3.77 




[4] 15940 


3.8 





librium internuclear separation R e , vibrational constant 
u! e ?» Yio and rotational constant B e w Yoi and those 
from theory. From [|[ the values for the state 3 0+ are 
taken. The values D e given in (5440 cm -1 ) and Q 
(5292 cm -1 ) are converted to T e by using the atomic level 
difference % — ^D 2 (20149.700 cm" 1 g|) and the D e 
value of the ground state (1081.64(2) cm" 1 [§|). 

The results for w e and B e from [6j show good agreement 
with the experimental results and deviate only by 2% and 
4%, respectively. 

A tabulated potential function (with high density of 
points) derived with the data from table IIII1 a highly 
precise potential describing only the unperturbed levels 
of this state in a point- wise and an analytic version (for 
definition see eq. @ below) can be found in the online 
material. 



B. The state FlL 

From all fluorescence progressions we derived in total 
153 levels of the state 1 1 II U . 

While fitting Dunham coefficients for the state l 1 ^ it 
soon became clear that this state is weakly perturbed. A 
small number of local perturbations is visible in the lad- 
der of f-levels obtained from the progressions of Q lines. 

Figure 2] shows the residuals for v' — 10 of the po- 
tential fit discussed later in this subsection and gives a 
typical picture of the observed perturbations. Each of 
the eight individual levels is derived from long fluores- 
cence progressions and the uncertainty for the energy is 
smaller than 0.01 cm -1 as also the Dunham fit in the 
last subsection indicated. Most of the f-levels shown are 
not significantly perturbed while there is one weak but 
significant local perturbation visible between J' = 158 
and J' = 182. The only e-level observed for this v' is 
more strongly shifted by a local perturbation caused by 
v' = 2 of the state 2 1 £+ as one can see using the data of 
section ITVCl 

Theoretically such perturbations of f-levels can only 
be caused by singlet states of type 1 II U , 1 A„ or 

triplet states with u-symmetry. Inspecting the ab initio 
calculations I6J the state 1 A U from the same asymptote 
X D + 1 S seems to be a good candidate as perturber, since 
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FIG. 4: (Color online) Fit residuals for v' — 10 of the state 
l 1 !!^ illustrating perturbations. 



TABLE V: Dunham coefficients Yik for the f-levels of the 
state ^Tlu and the reference isotopologue 88 Sr2. All values 
in cm' 1 . Yoo is calculated with respect to the potential min- 
imum of the ground state and includes the electronic term 
energy T e . 







16617.8249 0.02341446 -6.943x10" 

1 86.2991 -6.272xl0" 5 
9.7xl0~ 8 



2 -0.21086 

3 -2.3xl0~ 5 

4 -2.42xl0~ 5 



-3x10" 



it will have by far the largest Franck- Condon overlap with 
the state l 1 ^. 

Table IVl gives the resulting Dunham coefficients which 
can be used to reproduce the term energies of the unper- 
turbed levels from v' = to v' = 22 according to eq. ((T|). 
32 of the 153 observed levels, including all e levels, were 
removed from the Dunham fit because of perturbations. 
The weighted standard deviation of this fit is a = 0.85. 

Though all observed e-levels with rotational quantum 
number J' > 180 are systematically shifted to lower en- 
ergies than extrapolated from the observed f-levels, a 
consolidated description of the A-doubling can only be 
derived by a full deperturbation analysis of the system. 
The corresponding levels of the state 2 : E+ are shifted 
just into the opposite direction as the e-levels of the state 
l 1 ^. In the energy region, where for both states a suffi- 
ciently dense data set exists and precise frequencies can 
be calculated from the potentials, the extrapolated ener- 
gies of corresponding bands of both states reach the same 
values in the rotational ladder around J' = 250. The e- 
levels are additionally perturbed by the state 2 1 S+, but 
they could also be perturbed by highly excited vibra- 
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o.o 
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FIG. 5: (Color online) Residuals of the potential fit of the 
state plotted as a histogram. Few perturbed levels show- 
ing deviations of up to 1 cm -1 are out of scale in this figure. 



tional levels of the state 1 evaluated in the previous 
subsection. 

The Dunham coefficients were used to create a RKR 
potential which then served as a starting potential for the 
fit of an analytical potential representation, as defined in 
Q for the ground state. The central part of this potential 
(Ri < R < R a ) is represented as 



with the nonlinear mapping function 



(4) 



X = 



R — Rm 

R + bR„ 



(5) 



The inner repulsive wall (R < Ri) has the same form 
as for the potential used for the state l 1 !^ and is given 
in eq. ((2]) and similarly for the long range part (R > R a ) 
given in eq. ([3]). Here the a,j and T m are fit parameters 
while b, Ri, R m , R a and n are manually chosen to opti- 
mize the result. A and B from eq. ([2]) and the Ci from 
eq. © are adjusted with i = 8 and 10 to get continu- 
ously differentiable connections at the points Ri and R a , 
respectively. Uoo is calculated by adding to the difference 
of the atomic levels % - X B 2 of 20149.700 cm" 1 [H the 
ground state dissociation energy of 1081.64(2) cm -1 
The long range connection is mainly done to show the 
proper correlation of this state to the atomic asymptote. 

For the potential fit the same data set was used as 
for the Dunham fit. Figure [5] shows a histogram of the 
residuals of the potential fit of the levels assumed as be- 
ing unperturbed. The weighted standard deviation of 
the potential fit is a = 0.86. One can clearly see that 



TABLE VI: Potential coefficient for the state l 1 !!*,. The en- 
ergy is calculated with respect to the potential minimum of 
the ground state. The coefficients Cg, and Cio are simply 
used as continuously differentiable extensions from the deeply 
bound potential region to the atomic asymptote 1 So + 1 D2. 





-5.715X10" 1 cm" 1 


02 


1.6825387xl0 4 cm" 1 


a 3 


1.320905 xlO 4 cm" 1 




-1.03130xl0 3 cm" 1 


05 


-9.40946xl0 3 cm" 1 


<> a 


-1.92645xl0 4 cm" 1 


07 


-7.423611 xlO 4 cm" 1 


as 


-1.44756x10 s cm -1 


b 


-0.54 


Rm 


4.047308 A 


T m 


16617.8639 cm -1 


Ri 


3.524 A 


n 


6 


A 


1.4062979 xlO 4 cm" 1 


B 


8.341380xl0 6 cm" 1 !" 


Ra 


4.88 A 


C 8 


1.8205183xl0 9 cm" 1 A 8 


Cio 


-2.266162 xlO 10 cm" 1 A 10 


E/oo 


21231.34(2) cm -1 



the overwhelming amount of levels is centered around 
the middle with a full width at half maximum of less 
than 0.02 cm -1 = 2a and only few perturbed levels show 
larger deviations. 

In table lVII the resulting potential coefficients are listed 
which can be used to construct the potential according 
to eqs. (|![5]). 

A comparison of the potential from this work with the 
ab initio results is shown in figure |6l Here it is clearly 
visible that all calculations significantly overestimate the 
binding energy (i.e. too low T e values) of this state and 
again the earliest work Q shows by far the best agree- 
ment, a deviation of 375 cm" 1 , while the more recent 
works [4j, |7[ show deviations of about 800 cm -1 . Table 
IVIII gives a comparison of the spectroscopic constants 
showing significant deviations. Looking into the poten- 
tial plot given in Fig. 2b of [3] one can see that the adia- 
batic potential curve (state 6 l u ) has an avoided crossing 
with a ft = l u component of a triplet state ( 3 D + 1 S) 
at 5 A and thus the diabatic potential curve of the state 
l 1 n u to X D + X S would even deviate more strongly. The 
pronounced avoided crossing at the inner potential wall 
predicted by p and Q is also not observed, but if the 
binding energy of the state 2 1 n u by [f| is only little more 
overestimated than the binding energy of the state l 1 n„ 
it would be outside the investigated energy region. 

A point-wise potential of the state l 1 !!^ is given in the 
additional online material. 
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TABLE VII: Comparison of spectroscopic constants of the state 1 H u , the energy T e refers to the potential minimum of the 
ground state. 

Source T e [cm- 1 ] R e [A] w e [cm -1 ] B e [cm" 1 ] 
This work 16617.86(2) 4.0473(2) 86.300(3) 0.023415(2) 



[6] 


16243 


3.952 


96 


0.0246 


17] 


15821 


3.97 






1 


15818 


3.9 








3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 



R[A] 



FIG. 6: (Color online) Comparison of the potential energy 
curve with ab initio calculations of the lowest two 1 U n states 
from Boutassetta et.al. ml [Ic|. Czuchaj et.al. [7] (l 1 !^ and 
2 II„) and Kotochigova [J] (6(l u ) and 7(1^)). The energies 
are calculated with respect to the potential minimum of the 
ground state. The region with precise spectroscopic data is 
indicated. 

C. The state 2(A) J E+ 

For state 2(A) 1 £+ we collected the largest set of data, 
namely 1044 levels. Table I Villi gives Dunham coeffi- 
cients. They can be used according to eq. ([T]) to cal- 
culate the term energies of the unperturbed levels of all 
isotopologues up to v' = 23. The Dunham coefficients re- 
produce with a weighted standard deviation of a = 0.78 
the 1029 observed levels in the range up to v' = 23 except 
20 significantly perturbed ones. 

Compared to our earlier result Q where only data up 
to v' = 12 were available the potential function could 
be significantly extended to the range from 3.7 to 5.3 
A. Additionally to the 1009 levels used for the Dunham 
fit, 3 levels with v' = 28 and v' = 31 excited with the 
514 nm line of the Ar + laser are included in the poten- 
tial fit, of which one of the two levels with v' — 28 had 
to be removed again, likely because of a perturbation. 
The remaining 12 of the observed 1044 levels belong to 
higher vibrational levels and are not yet assigned to a v' 
quantum number. Histograms of the residuals of the fit 
are given in figure [7] for the three most abundant isotopo- 



30-, 




-0.04 -0.02 0.00 0.02 0.04 

Residuals [cm 1 ] 



FIG. 7: (Color online) Histograms of the residuals of the po- 
tential fit for the state 2(A) for the three most abundant 
isotopologues 86 Sr 88 Sr, 87 Sr 88 Sr and 88 Sr 2 . 

logues 86 Sr 88 Sr, 87 Sr 88 Sr and 88 Sr 2 . A few perturbed lev- 
els showing deviations up to 0.6 cm" 1 are not included. 
One can clearly see that the unperturbed or weakly per- 
turbed levels are centered around zero with a width at 
half maximum of less than ±0.01 cm" 1 which is the esti- 
mated measurement uncertainty for most levels and cor- 
relates to a weighted standard deviation of the fit of the 
1011 weakly perturbed levels of a — 0.88. The few levels 



10 



TABLE VIII: Dunham coefficients for the unperturbed levels up to v' = 23 of the state 2(A) 1 £+ for the reference isotopologue 
88 Sr2. All values in cm -1 . The energy refers to the minimum of the ground state potential. 

II k -> 1 2 3 

17358.7246 0.0219695 -5.98xl0" 9 

1 84.2157 -6.787xl0~ 5 -4x10"" -1.4xl0" 16 

2 -0.2665 -4.82xl0" 7 -1.7xl0" 12 

3 -0.00113 -7xl0" 9 

4 -l.lxlO" 5 
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FIG. 8: (Color online) The residuals of the 2(A) 1 E+ potential 
plotted in dependency of J' for v' — 13 and v' — 20 of the ref- 
erence isotopologue 88 Sr2 to illustrate perturbations between 
J" = 41 and 67 and J' > 250. 



of the weak isotopologues 86 Sr 2 ( 84 Sr 88 Sr) and 86 Sr 87 Sr 
show deviations from the predictions by the potential of 
less than 0.01 cm -1 , a quite satisfactory situation. 

Figure [5] illustrates observed perturbations by display- 
ing residuals as functions of J' for the vibrational levels 
v' = 13 and v' — 20 of the reference isotopologue 88 Sr2. 
For v' — 13 one can see a fast growing deviation for 
J' > 200 indicating a perturbation at very high J'. By 
comparing calculations for level energies of v' = 13 with 
those for the state l 1 !^ one finds that the v' = 21 stack 
of this state crosses the one of the state 2 1 £+ coming 
from smaller energies at J' = 271. The highest observed 
level of v' — 13 of 2 1 £+ corresponds to J' = 265. Such 
perturbations at high J' (moving slowly to lower J' for 
growing v') are observed for all v' of state 2 4 E+ for which 
a sufficient amount of data is available. For v' = 20 
one sees in figure [5] a perturbation between J' = 41 and 
J' = 67, which could stem from v' = 29 of state l 1 ^, 
but such extrapolation out of the known range of 1 n u 
is not yet reliable. 

In tabic ILXI the potential coefficients for the state 
2(A) 1 £+ are collected, which can be used to calculate 
the potential energy curve by using the eqs. (HI©. Here 



TABLE IX: Potential coefficients for the state 2(A) 1 £+. The 
energy is referenced to the potential minimum of the ground 
state. The long range coefficient C3 is the experimental value 
from [23l |. while the coefficients Cs and C10 are used for a 
continuously differentiable connection at the point R a . 



Ill 


3.89xl0~ 2 cm" 1 


111 


3.6222297xl0 4 cm -1 


az 


7.47647 xlO 2 cm" 1 


<i 4 


-5.9774017xl0 4 cm" 1 


as 


-1.2795502 xlO 5 cm" 1 


06 


-8.192268xl0 4 cm" 1 


a 7 


3.309500xl0 5 cm" 1 


as 


8.88141 xlO 4 cm" 1 


ag 


4.626819xl0 5 cm" 1 


b 


-0.33 


Rm 


4.17828223 A 


T 


17358.7496 cm" 1 


Hi 


3.49 A 


n 


6 


A 


1.4357908xl0 4 cm" 1 


B 


1.153771543 xlO 7 cm" 1 !" 


Ra 


6.1 A 


c 3 


5.9712(42) xlO 5 cm" 1 A 3 [23] 


Cs 


-2.7292900xl0 9 cm" 1 A 8 


C10 


5.5935384xl0 10 cm" 1 ! 10 




22780.122(20) cm" 1 



the long range coefficient C3 is taken from the photoas- 
sociation work |23;], while the coefficients Cs and C\o 
are simply derived for producing a continuously differ- 
entiable connection to the inner part of the potential at 
the point R a and to bridge the gap between this point 
(at 6.1 A) and the real long range region. Uoo is calcu- 
lated by adding to the frequenc y of the atomic transition 
% — ^Pi of 21698.482 cm" 1 [H] the ground state dis- 
sociation energy of 1081.64(2) cm" 1 [§]. 

Figure [5] shows a comparison of the resulting poten- 
tial to the ab initio calculations and table |X] summa- 
rizes the spectroscopic constants. For this electronic state 
the deviations are significantly smaller compared to the 
two other studied electronic states. Here the result from 
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TABLE X: Comparison of the spectroscopic constants of the state 2(A) 1 Ej with ab initio calculations [1, @ 5 0l an d an ear her 
experimental work The energy T e refers to the potential minimum of the ground state. 



Source 


T e [cm" 1 ] 


Re [A] 


aj e [cm : ] 


B e [cm" 1 ] 


This work 


17358.75(1) 


4.1783(1) 


84.215(1) 


0.021969(1) 


[15] 


17357.9(2) 


3.952 


85.07(34) 


0.02456(5) 


[6] 


17541 


4.099 


83 


0.0229 


\7] 


17290 


4.10 






1 


17269 


4.03 


88 
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FIG. 9: (Color online) Comparison of the potential for the 
state 2(A) with the ab initio calculations from Boutas- 
setta et.al. [f|, Czuchaj et.al. Q, and Kotochigova [J]. 

Czuchaj et.al. Q shows the best agreement and over- 
estimates the dissociation energy by only 69 cm -1 and 
the calculation by Boutassetta et.al. @ contrary to the 
other states gives the largest deviations by underestimat- 
ing the binding energy by 182 cm -1 . Table IXl indicates 
significant deviations of the present work to an earlier ex- 
perimental work [lol mainly for the rotational constant. 
They are due to a wrong rotational assignment which was 
already pointed out in our earlier publication [8j. 

A point-wise version of the potential and a list of the 
energy levels used for the potential fit can be found in 
the additional online material. 



V. CONCLUSION 

In this work we present a potential energy curve around 
the potential minimum for the state l 1 ^^ of which only 
little was known from matrix spectroscopy [lll - ll3l ] , a pre- 
cise potential for the minimum of the weakly perturbed 
state 1 II U) which was completely unknown experimen- 
tally, and a substantially improved potential energy curve 
compared to our earlier work Q for the state 2 1 £+. The 
results are compared to the currently available ab ini- 



tio calculations 0, S 0| which deviate strongly in several 
cases from our potentials and among each other. 

The binding energy of the state l 1 !!^ is not very well 
predicted by all available calculations, the deviations 
reach from an overestimation of the binding energy by 
about 400 cm -1 in Q to an underestimation by 3000 
cm -1 in 0, 0]- The presented potential is determined 
by experimental data in the range of v' = 1 to v' = 20 
which corresponds to internuclear distances from 3.5 A 
to 5.0 A and an energy interval of 1500 cm -1 . The po- 
tential function of the state l 1 shows a bent like an 
avoided crossing which corresponds to the similar case in 
Ca 2 @ where the coupling between and 3 IT U was 
studied. Thus the potential presented in table Mil should 
represent the lower adiabatic component — 0+ of the 
coupled system of these states. Surprisingly, all ab ini- 
tio calculations give the desired state 3 H U just around the 
energy where we find the bent in the potential curve. But 
Czuchaj et.al. 7] and Kotochigova [4] failed to predict 
the coupling between l 1 ^ and l 3 n u because the cross- 
ing of these states is high up in the repulsive branch. The 
calculation from Boutassetta et.al. Q put the crossing at 
the right ballpark, but the precision is far not sufficient 
to use their calculation for a start of a depcrturbation 
analysis. We will collect more data in the region where 
the strong coupling is expected for making a quantitative 
analysis of the coupled system l 1 ^^ + l 3 n u . 

The state l 1 !!^ is overestimated by all calculations in 
its binding energy, where @ shows the closest result de- 
viating by 375 cm -1 , while the recent works [1,0] deviate 
by almost 800 cm . Our experimental potential energy 
curve is determined by observed levels up to v' = 22 and 
covers the region of internuclear distances of 3.5 A to 
4.9 A and an energy up to 1800 cm -1 above the min- 
imum. A strongly pronounced avoided crossing at the 
inner repulsive potential wall with the state 2 1 IT M , pre- 
dicted by d, 0] within the observed energy region, is not 
observed. However, if the binding energy of the state 
2 1 n„ were overestimated a bit more than the binding 
energy of the state l 1 !^ in their calculations the avoided 
crossing would move outside the observation region. 

The electronic state best described by the theoretical 
calculations is the state 2(A) 1 £+, where all calculations 
find a binding energy deviating not more than 200 cm -1 
from our spectroscopic result. Q shows the best agree- 
ment by overestimating the binding energy by only 62 
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cm . However there are significant shifts in the inter- 
nuclear distances. The potential energy curve we present 
for this state is constructed from our data including vi- 
brational levels from v' = to v' = 31 in an interval of 
internuclear distances from 3.7 A to 5.3 A and up to 2400 
cm -1 above the potential minimum. 

In summary, for the currently available ab initio cal- 
culations d, H, 0] we found that surprisingly the earliest 
results @ are the most useful hints for our spectroscopic 
work and their potential curves were partly used during 
this work for first calculations of Franck-Condon factors 
and as a guide for choosing the excitation regions of a se- 
lected electronic state. Nevertheless, all calculations are 
far off from being applicable for direct calculations of ex- 
citation frequencies for experiments with cold molecules. 
As already mentioned for Ca2 [5j we found surprisingly 
large deviations in the electronic structure calculated in 
ab initio work for Sr2, too, and hope that our results 
might stimulate new efforts of the theorists. 

What can be concluded from our investigations regard- 
ing production of cold molecules from ensembles of cold 
atoms? By far the level v' = 4 of the state 2(A) 1 £+ has 
the best Franck-Condon overlap with ground state levels 
close to the asymptote and deeply bound ones. No other 
excited level of the investigated states is better suited 



than the one used in [9| for the investigation of the ground 
state asymptote. 

In P-Q it is proposed, based on the calculations in 
[H, to use levels of the state l 3 n 0n (1 0+) as interme- 
diate levels for a coherent transfer process of ultracold 
molecules from near asymptotic ground state levels to 
deeply bound levels of the ground state. At this moment 
we can already state that the potential energy curve of 
1 0+ looks completely different than the one calculated 
in [4( • If there are nevertheless rovibrational levels of the 
state 1 3 IT„ which are well suited for the proposed exper- 
iments we cannot judge at the moment since we did not 
reach that proposed frequency region. But we will try 
to answer this question when we continue on the coupled 
system of US^" - l 3 no u in the next future similarly to 
our study on Ca, 2 p|. 
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